Geometric phases for mixed states in interferometry 
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We provide a physical prescription based on interferometry for introducing the total phase of a 
mixed state undergoing unitary evolution, which has been an elusive concept in the past. We define 
the parallel transport condition that provides a connection-form for obtaining the geometric phase 
for mixed states. The expression for the geometric phase for mixed state reduces to well known 
formulas in the pure state case when a system undergoes noncyclic and unitary quantum evolution. 
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When a pure quantal state undergoes cyclic evolution 
the system returns to its original state but may acquire 
a phase factor of purely geometric origin. Though this 
was realized in the adiabatic context [Q, the nonadia- 
batic generalization was found in 0]. Based on Pan- 
charatnam's || earlier work, this concept was generalized 
to noncyclic evolutions of quantum systems Q. Subse- 
quently, the kinematic approach and gauge potential 
description ||[7j of geometric phases for noncyclic and 
non-Schrodinger evolutions were provided. The adiabatic 
Berry phase and Hannay angle for open paths were in- 
troduced ^ and discussed IjJ]. The noncyclic geometric 
phase has been generalized to non-Abelian cases Jl(J . Ap- 
plications of geometric phase have been found in molec- 
ular dynamics fll|| , response function of many-body sys- 
tem jl^,[l3| , and geometric quantum computation ]lj],[l5| . 
Noncyclic geometric phase for entangled states has also 
been studied In all these developments the geomet- 
ric phase has been discussed only for pure states. How- 
ever, in some applications, in particular geometric fault 
tolerant quantum computation |l4|,[l5|], we are primarily 
interested in mixed state cases. Uhlmann was probably 
the first to address the issue of mixed state holonomy, but 
as a purely mathematical problem In contrast, 

here we provide a new formalism of geometric phase for 
mixed states in the experimental context of quantum in- 
terferometry. 

The purpose of this Letter is to provide an opera- 
tionally well defined notion of phase for unitarily evolving 
mixed quantal states in interferometry, which has been 
an elusive concept in the past. This phase fulfills two 
central properties that makes it a natural generalization 
of the pure case: (i) it gives rise to a linear shift of the in- 
terference oscillations produced by a variable U(l) phase, 
and (ii) it reduces to the Pancharatnam connection |3| for 
pure states. We introduce the notion of parallel transport 
based on our defintion of total phase. We moreover intro- 
duce a concept of geometric phase for unitarily evolving 



mixed quantal states. This geometric phase reduces to 
the standard geometric phase for pure states under- 
going noncyclic unitary evolution. 

Mixed states, phases and interference: Consider a 
conventional Mach-Zehnder interferometer in which the 
beam-pair spans a two dimensional Hilbert space TL = 
{|0), |1)}. The state vectors |0) and |1) can be taken as 
wave packets that move in two given directions defined 
by the geometry of the interferometer. In this basis, we 
may represent mirrors, beam-splitters and relative U(l) 
phase shifts by the unitary operators 
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respectively. An input pure state /5; n = |0)(0| of the in 
terferometer transforms into the output state 

Pout = 

__ 1 / 1 + cos x 
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that yields the intensity along |0) as / cx 1 + cos%. Thus 
the relative U(l) phase x could be observed in the output 
signal of the interferometer. 

Now assume that the particles carry additional internal 
degrees of freedom, e.g., spin. This internal spin space 
Hi = C N is spanned by the vectors \k), k — 1,2, . . .N, 
chosen so that the associated density operator is initially 
diagonal 



Po 
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with u>fc the classical probability to find a member of 
the ensemble in the pure state |fc). The density operator 
could be made to change inside the interferometer 



po 



(4) 
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with Ui a unitary transformation acting only on the inter- 
nal degrees of freedom. Mirrors and beam-splitters are 
assumed to leave the internal state unchanged so that 
we may replace Um and Ub by Um = Um ® li and 
XJb = Ub ® li, respectively, li being the internal unit 
operator. Furthermore, we introduce the unitary trans- 
formation 
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The operators Um, Ub, and U act on the full Hilbert 
space H ® Hi- U correponds to the application of Ui 
along the |1) path and the U(l) phase \ similarly along 
|0). We shall use U to generalize the notion of phase to 
unitarily evolving mixed states. 

Let an incoming state given by the density matrix 
gin = pin® Po = |0)(0| <8 po be split coherently by a 
beam-splitter and recombine at a second beam-splitter 
after being reflected by two mirrors. Suppose that U 
is applied between the first beam-splitter and the mir- 
ror pair. The incoming state transforms into the output 
state 



g out = UBUjtfUUflftnUt, UtU^Uj, 
Inserting Eqs. (|l|) and (||) into Eq. (||) yields 
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The output intensity along |0) is 

/ cx Tr (u iPo U? +po + e- tx Uipo + e ix p U} 
cx 1 + |Tr {U iPo ) | cos [ X - argTr {U iPo )] , 



(8) 



where we have used Tr(poU^) = [Tr (Upo) 

The important observation from Eq. (g) is that the 
interference oscillations produced by the variable f (1) 
phase x * s shifted by <p — arg Tr(UiPo) f or any internal 
input state po, be it mixed or pure. This phase difference 
reduces for pure states po = \ipo)(ipo\ to the Pancharat- 
nam phase difference between Ui\ipo) and \ipo}- These two 
latter facts are the central properties for (f> being a natural 
generalization of the pure state phase. Moreover the vis- 
ibility of the interference pattern is v = |Tr (Upo) | > 0, 
which reduces to the expected v = | ('0o|£A|V'o) I for pure 
states. 

The output intensity in Eq. (^) may be understood as 
an incoherent weighted average of pure state interference 
profiles as follows. The state k gives rise to the interfer- 
ence profile 



I k oc 1 + v k cos [x - 4>k] ■ 



(9) 



where = |(fc|?7i|fc)| and 4>k 
the total output intensity 



&rg(k\Ui\k) . This yields 



I = ^ w k I k cx 1 + ^2 w k v k cos [x 



(10) 



which is the incoherent classical average of the above 
single-state interference profiles weighted by the corre- 
sponding probabilities w k - Eq. ( |lC)| ) may be written in 
the desired form 1 + v cos(x — <fi) by making the identifi- 
cations 
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Parallel transport condition and geometric phase: Con- 
sider a continuous unitary transformation of the mixed 
state given by p(t) — U (t) poW (t) . (From now on, we 
omit the subscript "i" of U.) We say that the state 
of the system p(t) acquires a phase with respect to pa 
if argTr[[/(t)/Oo] is nonvanishing. Now if we want to 
parallel transport a mixed state p(t) along an arbitrary 
path, then at each instant of time the state must be 
in-phase with the state at an infinitesimal time. The 
state at time t + dt is related to the state at time t as 
p{t + dt) = U(t + dt)U^(t)p(t)U(t)U^(t + dt). There- 
fore, the phase difference between pit) and pit + dt) is 
axgTr[p(t)U(t + dt)W(t)}. We can say p(t) and p(t + dt) 
are in phase if Tr[p(t)U(t + dt)U*(t)] is real and posi- 
tive. This condition can be regarded as a generalization 
of Pancharatnam's connection from pure to mixed states. 
However, from normalization and Hermiticity of p(t) it 
follows that Tr[p(t)U(t)U^(t)] is purely imaginary. Hence 
the above mixed state generalization of Pancharatnam's 
connection can be met only when 



Tr[p(t)U(t)U\t)] = 0. 



(12) 



This is the parallel transport condition for mixed states 
undergoing unitary evolution. On the space of den- 
sity matrices the above condition can be translated to 
Tr[p dU W] — 0, where d is the exterior derivative on 
the space of density operators. However, p(t) determines 
the N x N matrix U(t) (N being the dimension of the 
Hilbert space) up to N phase factors, and the single con- 
dition Eq. ( |l2|) while necessary is not sufficient to de- 
termine U(t). These N phase factors are fixed by the N 
parallel transport conditions 

(k(t)\U(t)U\t)\k(t)) = 0, k = 1,2,... N, (13) 

where the |fe(i))'s are orthonormal eigenstates of pit). 
These are sufficient to determine the parallel transport 
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operator U(t) if we are given a non-degenerate density 
matrix pit). 

The parallel transport condition for a mixed state pro- 
vides us a connection in the space of density operators 
which can be used to define the geometric phase. Thus 
a mixed state can acquire pure geometric phase if it un- 
dergoes parallel transport along an arbitrary curve. One 
can check that if we have a pure state density operator 
p(t) — %p(t))(tp(t)\ then the parallel transport condition 
Eq. (112) reduces to {ip(t)\tp(t)) — as has been discussed 
in 7|. p^ , p0{ which is both necessary and sufficient. 

Now we can define a geometric phase for mixed state 
evolution. Let the state trace out an open unitary curve 
r : t G [0,t] — > p(t) = U(t)p W(t) in the space of 
density operators with "end-points" p(0) — po and p(r), 
where U(t) satisfies Eq. (|l2|). The evolution need not 
be cyclic, i.e. p{r) ^ p . We can naturally assign a 
geometric phase 7 g [r] to this curve once we notice that 
the dynamical phase vanishes identically. The dynamical 
phase is the time integral of the average of Hamiltonian 
and can be defined as 



Id 



1 



dt Tr[p(t)H(t)] 



= -i i dt Tr[p {7t(t)f7(t)]. 
'o 



(14) 



Since the density matrix undergoes parallel trans- 
port evolution the dynamical phase vanishes identically. 
Moreover, the parallel transport operator U(t) should ful- 
fill the stronger condition Eq. (p"3|). Thus the geometric 
phase for a mixed state is defined as 



7s [r] =</> = argTr[ / 9 t/(t)]=arg 



w k v k e l<ih 

\ k 



(15) 



where exp(i/3fc) are geometric phase factors associated 
with the individual pure state paths in the given ensem- 
ble. The above geometric phase can be given a gauge 
potential description such that the line integral will give 
the open path geometric phase for mixed state evolution. 
Indeed the mixed state holonomy can be expressed as 



7fl [T]= / dt iTr[p W\t)W(t)} 



where 



= J iTi [poWUw] = J dfl, 
Tr[p a W(t)} 



W(t) 



|Tr[p C/t(t)]| 



U{t). 



(16) 



(17) 



and U(t) satisfies (|i]|. The quantity Q = iTi [p Q WUW] 
can be regarded as a gauge potential on the space of den- 
sity operators pertaining to the system. 

The geometric phase defined above is manifestly gauge 
invariant, does not depend explicitly on the dynamics 



but it depends only on the geometry of the open unitary 
path r in the space of density operators pertaining to the 
system. It is also independent of the rate at which the 
system is transported in the quantum state space. The 
geometric phase Eq. ( |l6| ) can also be expressed in terms 
of an average connection form 



7 S [r] 



/ S~] w k i(x k \dx k ) = / V^Wfefife, 
Jr l Jr t 



(18) 



where n k is connection-form and \xk(t)} — W(t)\k) is 
the "reference-section" for fcth component in the ensem- 
ble. To be sure, what we have defined is consistent with 
known results, we can check that this expression reduces 
to the standard geometric phase |^| fjj 

lg [Y] = arg<^(0)|V(r)) = f dt i(x(t)|*(*)) (19) 

Jo 

for a pure state p(t) — \ip(t)} (ip(t)\ when it satisfies par- 
allel transport condition. Here, \x(t)) 1S a reference state, 
which gives the generalised connection one-form (^J^). 

Purification: An alternative approach to the above re- 
sults is given by lifting the mixed state into a purified 
state \^>) by attaching an ancilla. We can imagine that 
any mixed state can be obtained by tracing out some de- 
grees of freedom of a larger system which was in a pure 
state 



l*> =Y,V^\k)s\k) ai 



(20) 



where \k) a is a basis in an auxiliary Hilbert space, de- 
scribing everything else apart from the spatial and the 
spin degrees of freedom. The existence of the above pu- 
rification requires that the dimensionality of the auxiliary 
Hilbert space is at least as large as that of the internal 
Hilbert space. If |\E') is the state at time t = and 
it is transformed to by a local unitary operator 

U(t) = U s {t) <g> I a then 



!*(*)> =^2V^kU s (t)\k) s \k) a . 



(21) 



The inner-product of initial and final state 



(*(0)|tf(t)) =J2^(k\U(t)\k) = Tr(C/(t)p ) (22) 

k 

gives the full description of the modified interference. 
Indeed by comparing Eqs. (||) and (|22]), we sec that 
arg(#(0)|#(t)) is the phase shift and |(#(0)|#(i))| is the 
visibility of the output intensity obtained in an interfer- 
ometer. 

The parallel transport condition, given by Eq. ( |l2| ) , fol- 
lows immediately from the pure state case when applied 
to any purification of po, i.e. 
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= {#(t)\*(t)) = ^w k {k\Ut(t)U(t)\k) 
k 

= Tr\porf(t)U(t)] = Tr[p(t)U(t)uHt)]. (23) 

Thus a parallel transport of a density operator p(t) 
amounts to a parallel transport of any of its purifications. 

Example: Consider a qubit (a spin-i particle) whose 
density matrix can be written as 



P : 



(24) 



where f is a unit vector and r is constant for unitary evo- 
lution. The pure states r — 1 define the unit Bloch sphere 
containing the mixed states r < 1. Suppose that dur- 
ing the time evolution f traces out a curve on the Bloch 
sphere that subtends a geodesically closed solid angle Ct 
Ipjf . The two pure states |±; f • cr) acquire noncyclic geo- 
metric phase =pJ7/2 and identical visibility Vj r =V-=r\. 
Using Eq. ( |l5| ) we obtain the geometric phase for T 



arctan r tan 







= 7 9 [r] 



The visibility v = |Tr (Upo) | is given by 



v = rj\j cos^ — + r 2 sin — . 



(25) 



(26) 



For cyclic evolution we have rj — 1 but the mixed state 
v 7^ 1 due to the square root factor on the right-hand side 
of Eq. (|2rj|). Moreover Eqs. (25) and ( p6|) reduce to the 
usual expressions for pure states <fi — — f2/2 and v = rj 
by letting r = 1. 

In the case of maximally mixed states r = we obtain 
= argcos(f2/2) and v = |cos(fi/2)|. Thus the output 
intensity for such states is 



Joe 1 



= I 



cos — cos 
1 2 1 

CI 

cos — cos%. 



X — arg cos ■ 



(27) 



Early experiments |23|-|25|| to test the An symmetry of 
spinors utilized unpolarized neutrons. Eq. show that 
in these experiments the sign change for £1 = 2tt is a con- 
sequence of the phase shift <f> = arg cos 7r = tt. 

Note that 7 ff [r] in Eq. (|2^) equals the geodesically 
closed solid angle on the Poincare sphere iff r = 1 . In the 
mixed state case the geometric phase factor is weighted 
average of the solid angles subtended by the two pure 
state paths on the Bloch sphere. 

In conclusion, we have provided a physical prescription 
based on interferometry for introducing a concept of total 
phase for mixed states undergoing unitary evolution. We 
have provided the necessary and sufficient condition for 
parallel transport of a mixed state and introduced a con- 
cept of geometric phase for mixed states when it under- 
goes parallel transport. This reduces to known formulas 



for pure state case when the system follows a noncyclic 
and unitary quantum evolutions. We have also provided 
a gauge potential for noncyclic evolutions of mixed states 
whose line integral gives the geometric phase. We hope 
this will lead to experimental test of geometric phases for 
mixed states and further generalization of it to nonuni- 
tary and nonlinear evolutions. 
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